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Losses in superconducting planar resonators are presently assumed to predominantly arise from surface- 
oxide dissipation, due to experimental losses varying with choice of materials. We model and simulate the 
magnitude of the loss from interface surfaces in the resonator, and investigate the dependence on power, 
resonator geometry, and dimensions. Surprisingly, the dominant surface loss is found to arise from the metal- 
substrate and substrate-air interfaces. This result will be useful in guiding device optimization, even with 
conventional materials. 



Superconducting coplanar waveguide (CPW) res- 
onators are critical elements in photon detectioni, quan- 
tum computatioii^ri^, and creating and decohering quan- 
tum photon states^ii. Such applications are limited by 
the energy decay time. One prominent source of de- 
coherence at low powers in resonators has previously 
been found to be two-level states (TLSs) on the various 
surfaces^"^^. Knowing the TLS locations is important 
for improving the resonators. Previous measurements 
and simulations suggested that the exposed metal surface 
(metal-air interface) is a crucial decoherence source^rJi, 
which has driven research in using non-oxidizing super- 
conductors for quantum devices. However, we show with 
simulations and a model that, for typical metal oxide pa- 
rameters, the more likely source of CPW loss is instead 
the metal-substrate and substrate-air interfaces, thereby 
changing the approach needed to reduce losses. 

Thin-film resonators have three types of amorphous 
interfaces that can contain TLSs and thus introduce sur- 
face losses (Fig.[l]): metal-air (ma), metal-substrate (ms), 
and substrate-air (sa). Each interface can include a thin 
oxide or contaminant layer sandwiched between the two 




FIG. 1. (Color online) Coplanar waveguide dimensions and 
interfaces. Dimensions are the center width w, gap width g, 
metal height h, etch depth d, and assumed interface thickness 
i=3nm. This structure has metal-air (blue), metal-substrate 
(red), and substrate-air (black) interfaces. The 3nmx3nm 
corner square (yellow) is treated separately. 



primary layers. The total loss tangent for these thin lay- 
ers is X^jPi tan (Si, where surface interface type i has loss 
tangent tan 6i and participation ratio^ii^ 



P, 



= W~H,e, ds IE 



(1) 



for thickness ti, dielectric constant e^, length coordinate 
s, and energy per unit length W. Using the boundary 
conditions on the electric displacement gives the three 
interface participation ratios^^: 



ds IE,, 



P-maW/t-ma. = £„ 
P-msW/t^s = (Cs/ems) / ds \Esl_ 
PsaW/Us, = esa / ds\E. 



2 + e-i 



(2) 
(3) 

rfs|^a±P, (4) 
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where E^ {E^) is the electric field in the air (substrate) 
outside the interface and E\\ {E±) is the electric field 
component parallel (perpendicular) to the interface. 

Here, we take all dielectric constants to be of order 
e ^ 10; typical of metal oxides. Then, pma and Psa,.± are 
of order 1% of pms and Psa.jj- Thus, if all interface loss 
tangents and thicknesses are similar, the substrate-air 
and metal-substrate interfaces are 100 times more lossy 
than the metal-air interface. 

To accurately compare to this model, the participa- 
tion ratios were numerically calculated using the finite- 
element solver COMSOLl^ with Eqs. ^-^ used to ex- 
tract the participation ratios from surface fields^''. As 
shown in Table U pma is 40-60 times smaller than p,ns 
and psa, which corresponds to the difference in prefac- 
tors in Eqs. (I2])-(|4]), validating the discussion above. In 
addition, for typical interface parameters (t ~ 3 nro^, 
tan (J ^ 0.002^^), the metal-air interface gives a quality 
factor Q ^ 10^, much greater than typically measured 
values of 10^; the measured values are more similar to the 
values of 1.5 — 2.5 x 10^ predicted for the metal-substrate 
and substrate-air interfaces. 



TABLE I. Simulated losses for three resonator geometries^*. 
Losses are calculated for ts = 10 and surface dielectrics with 
e — 10, i=3nm, and loss tangent 0.002; dimensions are 
as in Fig.[T] Both coplanar waveguide (CPW) geometries, 
with identical dimensions, have similar metal-air (ma), metal- 
substrate (ms), and substrate-air (sa) losses, but etching the 
exposed substrate (Etched CPW) substantially reduces cor- 
ner (c) loss. A microstrip geometry (with dielectric height s) 
has significantly less ma and sa losses than CPWs. 



♦♦♦ma •••c ■■■ms •••sa 



Type 


Dimensions (/im) 


Loss xlO*^ 


w 


h 


9 


d 


ma 


ms 


sa 


c 


CPW 


5 


0.1 


2 





0.10 


6.13 


4.02 


1.32 


Etched CPW 


5 


0.1 


2 


0.01 


0.11 


4.64 


5.25 


0.39 


Microstrip 


20 


0.2 


3=2 





0.02 


6.60 


0.79 


0.38 



Participation ratios have dominate contributions from 
the edges around the coplanar gap^^; the scaling of the 
electric fields with distance is displayed in Fig.[2l For 
distances r from corners much greater than other relevant 
dimensions, the field scales as r~^/'^, the predicted scaling 
for the field from a flat edge^^. The electric fields on the 
metal-air interface and on the unetched metal-substrate 
interface scale as r~^/^, the predicted behavior for the 
field from a metal 90° corner—; in this case, 20% of the 
participation ratio is within 1% of the length at the edge. 
Because of this, we also calculate the participation ratio 
Pc for a 3 nm square at metal-air-substrate corners (as 
shown in Fig. [J), and find the loss from this small corner 
is only 3-4 times smaller than from the metal-air and 
substrate-air interfaces. 

Etching the substrate in the coplanar gap flattens the 
field dependence on r at the metal-air-substrate corner 
(Fig. [2]). This implies etching the substrate significantly 
reduces pc while leaving p,„a and p,„s unchanged, along 
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FIG. 2. (Color online) Electric field scaling on distances for 
CPW resonators, (a) is for an ordinary CPW geometry, while 
(b) is for a CPW with the substrate etched in the coplanar 
gap. Arrows indicate for which distance r from the corner the 
electric field has the indicated scaling. As expected^-, the field 
far from metal edges (b) scales as r~^ " and fields at the metal 
edge in both cases scale as r~^' ^; a substantial fraction of the 
energy in both cases is concentrated at the corners. Etching 
the substrate in the coplanar gap reduces this dependence 
over a length corresponding to the etch depth, reducing the 
corner participation ratio. 
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FIG. 3. (Color online) Geometric dependence of participation 
ratios for the metal-air (ma), corner (c), metal-substrate (ms), 
and substrate-air (sa) interfaces. We assume 3 nm surface di- 
electrics with Ema = Ems = Esa = 10 aloug with /i=100 um 
and d=10nm. (a) Open symbols are for g = w, where the 
loss decreases as l/w. Filled symbols are for fixed g =20 /im 
(indicated by dashed line), where the loss decreases as w~^'^. 
(b) Plot of loss vs. w for fixed w + 2g=40 fim. Minimum 
loss occurs at a characteristic impedance Zq =50-60 Q.. In 
all three cases, the participation ratios scale together, imply- 
ing that changing geometrical parameters can not determine 
which interface is dominant. 



with a potential decrease in psa due to lower surface fields. 
Remarkably, even a 10 nm etch reduces pc by 70%, while 
a 2 ^m etch reduces pc by 99% along with psa by 50%. 

The microstrip geometry also changes the participa- 
tion ratios, as shown in Table HI Compared to CPW 
resonators, this approach significantly reduces pc, Pma, 
and psa while leaving p^s unchanged (Table |I|; hence, 
microstrip resonators are especially useful if the metal- 
substrate interface is the least lossy interface. This differ- 
ence between the interfaces implies that it is possible to 
determine if the metal-substrate interface is dominant by 
comparing losses from CPW and microstrip resonators. 

One might expect varying CPW resonator dimensions, 
as shown in Fig. [31 would determine dominant interfaces. 
For a fixed coplanar gap g, the loss decreases as w^^'^ 
for w < g and flattens off for w > g. Since larger widths 
increase loss by producing slotline modes, increasing radi- 
ation, and allowing trapped flux in the center stripi^, we 
characterized the case w = g, where loss is proportional 
to l/w. If the distance w + 2g between the ground planes 
is kept fixed, then minimal loss occurs at Zq = 50 — 60 il. 
However, in all three cases, the participation ratios for 
all four interface types have nearly the same dependence. 
Hence, varying dimensions can reduce loss but makes de- 
termining the dominant interface difficult. 

One potential way to determine the key interfaces 
is by measuring the power dependence of the lossiS, 
since different interfaces have different field dependences 
from their corners (Fig. [2). For a 90° corner, where 
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FIG. 4. (Color online) Plot of surface participation ratio 
vs. electric field for different field dependences. The case 
E ~ r~^'^ , for a metal corner [Eq. ([5]), blue], has similar 
dependences if Tc = (solid) or if rc/ro — 10~* (dashed); 
both are similar to the TLS saturation theory case [E = Eq, 
black] and the simplified formula [Eq. ((7)l for a = 1.5, green]. 
For the case E r^ )"~^ , for a metal edge ]Eq. ([6], red), the 
Ta/rQ = 10~* curve (dashed) has a broader crossover, while 
the Tc = (solid) case logarithmically diverges. 



E = Eo{r/ro) ^^^, the surface participation ratio isH 



te 



3E^ro 




(5) 



where Es is the saturation field, ro is a characteristic 
length, and re is a lower cutoff. For a metal edge, where 
E = £'o(r/ro)^^/^, the surface participation ratio is^^ 



1 



te 



2E^ro log 



1 



E^ 
§ 



ro ^ S2 



(6) 



which gives a logarithmic divergence. Comparing these 
results in Fig.|4l we find the r~^/^ edge model has a much 
broader crossover with the the drive field, while the r~^/'^ 
corner model is much closer to the sharp crossover of the 
simple TLS theory^. These models are also well described 
(Fig.U by the experimentally-based fitting formulai^ 



P_ 
te 



cE^ro 



[1 + 0.9c{Eq/E,) 



q1 1/q ' 



(7) 



where c = 3[1 — (rc/roY^^]- Here, a ~ 1.5 for the r~^/^ 
case, similar to an experimentally-determined relation^, 
while a ~ 0.75 for the r~^/^ case; this indicates that a 
corner, not an edge, dominated the experimental loss. 

The surface loss model is partially validated by qualita- 
tively explaining experimental results observed by other 
groups. For instance, HF-terminating a silicon substrate 
before growing an Nb film reduced the TLS loss by a fac- 
tor of three compared to in situ RF cleaning, confirming 
the importance of the metal-substrate interface^^. While 
the roughness in the coplanar gap did not affect their TLS 



loss, they substantially overetched into the substrate, 
so the substrate-air roughness experienced substantially 
lower fields. In addition, it was found that etching into 
the substrate reduced loss^'^, confirming the importance 
of the substrate-air interface. 

While previous experiments suggested that the metal- 
air interface is dominant, these results can be reexplained 
in terms of the model presented here. Our group pre- 
viously found MBE-grown Re resonators had lower loss 
than sputter-grown Al resonators^, which was attributed 
to the metal-air interface. This may instead be from dif- 
ferences in surface preparation and cleanliness between 
the growth conditions, affecting the metal-substrate in- 
terface. In addition, Nb resonators on Si and sapphire 
(AI2O3) substrates had similar losses and were more lossy 
than Al, Re, or TiN resonators; this was claimed to be 
from oxide formation at the metal-air interfaceii. How- 
ever, their Nb-Al203 resonators were grown on A-plane 
AI2O3, with which oxygen atoms diffuse into the Nb 
filmic, making the metal-substrate interface significantly 
thicker and thus more lossy. In fact, others showed that 
Nb resonators on R-plane AI2O3 have losses similar to 
Re and Al^*', indicating that the loss was not from oxide 
at the metal-air interface. 

This model is not even restricted to CPW resonators. 
The Nb cavity resonators of Raimond, Haroche, and 
Brune'^^ have a loss of 1/Q = 2.4 x 10~-^^. For typ- 
ical parameters, the surface loss model gives a loss ^ 
tan (5,naPma/ema = 2.2 X 10^^^, closc to the experimental 
data. 

In conclusion, we have developed a model for the res- 
onator loss from interfaces. We find that, for realistic val- 
ues, the metal-substrate and substrate-air interfaces are 
dominant, with participation ratios of order 100 times 
that of the metal-air interface. The loss can therefore be 
reduced by improving the metal-substrate and substrate- 
air interfaces, using microstrips with clean dielectrics, 
and increasing dimensions. 

This work was supported by lARPA under ARO award 
W911NF-09-1-0375. M.M. acknowledges support from 
an Flings Postdoctoral Fellowship. 
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Calculations are provided for the equations in the manuscript "Surface loss simulations of superconducting 
coplanar waveguide resonators". We provide a table of surface loss participation ratios for different geometries. 
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I. DERIVATION OF SURFACE LOSS MODEL 

In Eq. (1) of the main manuscript, the participation 
ratio for interface i is given bj^*^ 



Pi 



W-'Ue, / ds \E 



(SI) 



where the interface has a small thickness U, dielectric 
constant e^, and length coordinate s and where the res- 
onator structure has an energy per unit length W. 

The metal-air (ma) interface consists of the metal, 
a thin metal oxide with thickness tma — 3 nm and di- 
electric constant ema, and the outer air (vacuum) with 
Ea = 1- The electric field must be perpendicular to the 
metal surface, and because the interface layer is thin, 
we also approximate it as perpendicular in the dielectric, 
so £'nia = -EiiiaX- The Continuity of eE at the metal- 
oxide and air interface requires ema£'maj_.t = Eg^±^f Since 
the oxide is thin, E does not change significantly over 
the oxide thickness. Combining all these results gives 
Ema. ~ Eg.^/cma, SO the participation ratio of the metal- 
air oxide is 



PmaW^/^ma = Ema / ds \Ema\ 

= Ema / ds |£^ai/emaP 

'"' hs\E,i_\\ (S2) 



For the metal-substrate interface, we assume a thin 
dielectric layer of unknown origin between the metal and 
substrate, which might arise from a chemical reaction 
of the metal to the substrate or chemi- or physi-sorbed 
water on the wafer surface. As before, the electric field 
is perpendicular to the metal and the continuity of the 
displacement field requires ems^ms±,t = Es±^t, where ms 
represents this dielectric and s the substrate. Thus, we 
find Ems ~ Es_\_es/ems, so the participation ratio of the 



metal-substrate layer is 

PmsW/tms = Cms / ds \Ems\'^ 

= ems / ds \Es±es/ems\' 



(S3) 



For the substrate-air interface, there can be a dielec- 
tric layer from surface water or other contaminants from 
the air. described by a subscript sa. In addition to the 
perpendicular electric field which obeys Esa ~ ^a±/esa 
as before, there are also parallel field components obey- 
ing the boundary condition E^^u = E^aii = E^u, since the 
interface layer is thin. Hence, the participation ratio of 
the substrate-air interface layer is 



PsaW/U, = e,a f ds (|S,a|||' + \Esa±\^) 

= e,,jds |Sa||P+e-a'/rfs \Ea 



(S4) 
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II. SIMULATION APPROACH 

The coplanar and microstrip structures were simulated 
using the electric quasi-statics component of the finite 
element solver COMSOL's AC/DC module'''. We sim- 
ulated a two dimensional cross-section with half of the 
resonator, using symmetry to account for the other half. 
We used adaptive meshing as a starting point and then 
performed additional meshing around the edges and the 
corners. 

To determine the participation ratios, we initially 
treated the interfaces as 3 nm thick dielectrics with di- 
electric constant e = 10. However, this area approach is 
computationally expensive since it requires meshing on 
the nanometer scale over distances of hundreds of mi- 
crons. As such, we primarily calculated the participa- 
tion ratios by computing the electric field on all bound- 
ary interfaces with the interface dielectrics excluded from 
the model and then applying Eqs. (|S2P - (jS4[ ). This is less 
computationally expensive as there is no thin dielectric 
layer explicitly included at the interfaces which needs to 
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TABLE SI. Simulation results for a variety of microwave resonators, obtained from the primary surface-based model [s], the 
area model [a], and a previous calculation*. Coplanar waveguide and microstrip resonator dimensions are as indictated in 
Fig. lS2l the angle 6 enclosed by the metal at the metal-substrate-air corner is assumed to be 90° unless otherwise noted. The 
+ sign in the metal-air column data is for a 3nm by 3nm area (c) at the intersection of the metal-air (ma), metal-substrate 
(ms), and substrate-air (sa) interfaces, and represents an entry that could be split among the three interface types. It is placed 
in the metal-air column since there it gives the greatest proportional uncertainty and because pma + Pc and not Pma alone is 
comparable between the area and surface models. We assume eg = 10 and surface dielectrics with e = 10, thickness 3nm, and 
loss tangent 0.002. 



type 


dimensions 


capacitance 
pF/m 


metal- air 
Pma (ppm) 


metal-sub. 
Puis (ppm) 


sub. -air 
Psa (ppm) 


loss metal-air 
xlO^ 


loss metal-sub. 
xlO** 


loss sub. -air 
xlO^ 


coplanar 


w,h,g,d 
















cl[a] 


5, 0.1, 2, 0.01 


162 


119+167 


2200 


2541 


0.24+0.33 


4.40 


5.08 


c2[s] 


5, 0.1, 2, 0.01 


162 


56+196 


2322 


2624 


0.11+0.39 


4.64 


5.25 


c3[a] 


5, 0.1, 2, 


163 


290+387 


2234 


2286 


0.58+0.77 


4.47 


4.57 


c4[s] 


5, 0.1, 2, 


163 


52+662 


3065 


2011 


0.10+1.32 


6.13 


4.02 


c5 [a*] 


5, 0.1, 2, 




600 




2000 


1.2 




4.0 


c6 [s] 


5, 0.1, 2, 2 


104 


44+6 


2690 


1032 


0.09+0.01 


5.38 


2.06 


c7[s] 


5, 0.025, 2, 0.01 


161 


55+209 


2376 


2735 


0.11+0.42 


4.75 


5.47 


c8[s] 


2, 0.1, 20, 0.01 


68 


33+111 


1394 


1594 


0.07+0.22 


2.79 


3.19 


c9 [s] 


5, 0.1, 20, 0.01 


85 


18+60 


847 


928 


0.04+0.12 


1.69 


1.85 


clO [s] 


5, 0.1, 20, 


85 


17+207 


1091 


764 


0.03+0.41 


2.18 


1.53 


ell [s] 


5, 0.1, 20, 0,61 = 45° 


169 


32+1414 


3727 


2267 


0.06+2.83 


7.45 


4.53 


cl2 [s] 


5, 0.1, 20, 0,61 = 135° 


158 


104+695 


2841 


1963 


0.21+1.39 


5.68 


3.93 


microstrip 


w, h, s, d 
















ml [s] 


20, 0.2, 2, 0.01 


985 


10+45 


3155 


526 


0.02+0.09 


6.31 


1.05 


m2 [s] 


20, 0.2, 0.2, 0.01 


8964 


7+55 


29942 


409 


0.01+0.11 


59.9 


0.82 


m3 [s] 


10, 0.2, 2, 0.01 


539 


19+82 


3301 


964 


0.04+0.16 


6.60 


1.93 


m4 [s] 


20, 0.02, 2, 0.01 


983 


10+49 


3185 


557 


0.02+0.10 


6.37 


1.11 


m5 [s] 


20, 0.2, 2, 


987 


9.3+189 


3301 


397 


0.02+0.38 


6.60 


0.79 


m6 [s] 


20, 0.2, 2, 2 


914 


4.6+1.9 


2924 


291 


0.009+0.004 


5.85 


0.58 


m7 [s] 


20, 0.2, 2, -2 


1006 


1.5+3.2 


3192 


241 


0.003+0.006 


6.38 


0.48 



be carefully meshed. As indicated in Fig.[ST] and Table 



for two different pairs of simulations, Pn 



Ps. 



and 



Pma + Pc as calculated by these two approaches typically 
agree to within 15%, although pma alone differs by a fac- 
tor of at least two. This means that the total metal-air 
interface includes an indeterminate significant fraction of 
the corners. 

We also assumed all surfaces were smooth. Simulations 
indicate that incorporating smooth bumps on the order 
of the interface thickness increase the participation ratios 
and thus loss by a factor of order unity. The value of 
this factor depends on the interface thickness and on the 
defect density. 



III. RESULTS FOR DIFFERENT GEOMETRIES 

Numerical results for a variety of coplanar and mi- 
crostrip resonator geometries are presented in Table ISII 
We have calculated the participation ratio pi and loss 
Pi tan (5 for a dielectric with thickness 3nm, dielectric 



constant e — 10, and loss tangent tan (5 = 0.002, typical 
values for metal or silicon oxides^. Since the participa- 
tion ratio is proportional to thickness, these values can 
easily be scaled for other parameters. Assuming these 
parameters, the loss from the metal-air interfaces is typi- 
cally below 10~^. The second quantity in the sum for the 
metal-air columns arises from the 3 nm by 3 nm corner at 
the metal/substrate/air interface. As this is a small area, 
it shows the sensitivity of the loss to this inside corner 
and indicates the uncertainty in the metal-air prediction. 

For coplanar resonators, a significant effect on total 
loss comes from etching into the substrate within the 
coplanar gap. This is important because the etching re- 
duces the divergence of the fields at the corner, as shown 
in Fig. 2 of the main paper, and because the etching re- 
duces the fields parallel to the metal-substrate interface 
due to the increased distance from the metal traces. Be- 
tween the pair of cases cl and c2 and the pair c3 and 
c4, it is apparent that pma is reduced by a factor of 2- 
3. Further deep etching (2/im, case c6) reduces pma by 
an additional factor of 5 and psa by a factor of 2 while 



S3 



leaving unchanged Pms- This change has been experimen- 
tally tested for Si substrates^, where a feature in the loss 
versus power saturation curve was identified as substrate- 
air loss. After etching the substrate caused the feature 
to disappear and resulted in half the loss, which is con- 
sistent with the halving of psa. in case c6. This result 
is consistent with the discussion in the main paper indi- 
cating that the substrate-air interface is a dominant loss 
mechanism. 

As experimentally seen by the Delft group^, different 
geometries of coplanar resonators result in somewhat dif- 
ferent losses. In comparing cases cl and c2, there is min- 
imal difference between the area and surface models ex- 
cept for some change in pma from the inside corner, and 
comparing cases c4 and c5 indicate that these data are 
similar to previous work from our group**. The metal 
thickness is shown to have little effect in case c7. Case c8 
shows that decreasing the width makes the loss increase. 
However, in case c9, increasing the gap from 2 to 20 /im 
gave roughly a factor of 3 reduction in all losses. 

Sloped sidewalls are also seen to give different losses by 
comparing cases cf 1, c4, and cl2, where the metal angle 
6 at the substrate-corner corner fFig. lS2[) was varied. All 
interfaces except the metal-air interface had the greatest 
loss in case ell, where the sidewall slope was 9 — 45°, and 
the least loss in case cl2, where the overetched sidewalls 
had 9 = 135°. This is expected since the electric field 
is predicted to scale with the distance r from corners 
as r-3/7 for 9 = 45°, r'^/^ for 9 = 90°, and r'^^^ for 
9 = 135°*', thus giving the least field divergence, and thus 
the lowest loss, at the overetched corner. The metal-air 
interface exhibits the opposite trend, which is consistent 
with the same argument for the top corner. 

Microstrip resonators show significantly higher capaci- 
tance per length, which for the same interface energy re- 
sults in lower loss. In the base case ml, psa and Pma +Pc 
are both much less than the corresponding values for 
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FIG. S2. Coplanar and microstrip dimensions, (a) Dimen- 
sions for coplanar waveguide simulations. We assume 6 = 90° 
unless otherwise specified, (b) Dimensions for microstrip sim- 
ulations. 



coplanar resonators. However, microstrip resonators also 
have a larger pms contribution, approximately equal to 
the distance ratio 2t/s for oxide thickness t. In fact, the 
thin dielectric of case m2 compared to the base case ml 
has a capacitance and Pms approximately 10 times that 
in case ml and similar psa and Pma- Another geometric 
parameter that is important is the width w, with which 
Pma and Psa scalc inversely, as shown by case m3. How- 
ever, as indicated in case m4, changing the metal height 
has minimal effect on the participation ratios. 

As with coplanar resonators, changing the depth of 
etching of the exposed substrate also has a significant 
effect on loss. The lack of dielectric etching in case m5 
results in an increase in pma from the inside corner of the 
metal, which is mostly compensated for by a decrease 
in Psa- Case m6 shows that a deeper etch significantly 
reduces both of these terms. If, instead of etching the 
exposed dielectric, the dielectric extends up the sidewall 
of the metal (as in case m7), the losses are similar to that 
of a deep etch. Overall, the improvements with etching 
the exposed dielectric for microstrip resonators are pri- 
marily at the metal-air interface, but some effect (factor 
of 2) is seen for the substrate-air interface. However, the 
most important concern for a microstrip geometry is a 
low loss metal-substrate interface, as this loss was domi- 
nant in all cases. 
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FIG. SI. Etched coplanar waveguide participation ratios. 
Participation ratios are given in parts per million (ppm) for 
the metal-air (ma), metal-substrate (ms), and substrate-air 
(sa) interfaces, as calculated with the surface fields approach 
with the metal-air-substrate corner (c) treated separately. 
Participation ratios for the area approach are shown sepa- 
rately. We assumed 3 nm surface dielectrics with ema = tms = 
Esa = 10; the geometry is given in Table [STI in cases cl and c2. 
The participation ratios Pma+Pc, Pms, and psa agree to within 
15%, validating the simulations against numerical errors. 



IV. DERIVATION OF POWER DEPENDENCE OF 
PARTICIPATION RATIOS 

One potential way to determine which interface dom- 
inates the loss is by measuring the power dependence of 
the losS". When the saturation of surface two-level states 
(TLSs) at the field Eg is considered, the surface partici- 
pation ratio p is given by 



P_ 
te 



E^dr 



y/T+WjE^ 



(S5) 



where the interface has thickness t and dielectric constant 
e. The electric field is assumed to be dominated by a 
feature such as a corner with length coordinate r from 
this feature, characteristic length tq, and cutoff length 
re < tq. 

For a square corner, the field scales® as E = 
Eo{r/ro)~'^^^. so substituting this into Eg. (jS5|) gives a 
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surface participation ratio of 



P _p2 



= SE^ro 



(ro/r)^/^ dr 



^/i + MMW^hW 




(S6) 



For a thin edge at distances much greater than the film 
thickness, the field scales^ as £' = £'o(r/ro)~^/^, so sub- 
stituting this into Eq. (|S5ll gives a surface participation 
ratio of 



ra 



(ro/r) dr 



te °7.^ ^l + iE'^/Emro/r) 



2Elro log 



El 



re I I re j_ EsL 

7^ + V7;r + 4 



(S7) 
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